APS/123-QED 



Point processes and stochastic displacement fields 
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The effect of a stochastic displacement field on a statistically independent point process is ana- 
lyzed. Stochastic displacement fields can be divided into two large classes: spatially correlated and 
uncorrelated. For both cases exact transformation equations for the two-point correlation function 
and the power spectrum of the point process are found, and a detailed study of them with important 
paradigmatic examples is done. The results are general and in any dimension. A particular attention 
is devoted to the kind of large scale correlations that can be introduced by the displacement field, 
and to the realizability of arbitrary "superhomogeneous" point processes. 
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I. INTRODUCTION 

Point processes (i.e., stochastic spatial distributions of 
point-particles with identical mass) are very useful math- 
ematical models of many n— body and complex systems. 
Crystals (regular, perturbed, and/or defected) 0,0,0, 
quasi-crystals structural glasses, fluids 3, cosmolog- 
ical self-gravitating systems [a, , and also computer im- 
age processing problems @, and bio- metrical studies 
are only some examples of systems which are usually rep- 
resented as specific point processes with appropriate spa- 
tial correlation properties. 

The study of this branch of stochastic processes and 
the discovery of new statistical properties can be of fun- 
damental importance in many scientific topics. Many 
mathematical studies have been already done about this 
class of processes and many useful results have been de- 
rived (e.g., see [13 

One important question about a point process is what 
happens to its statistical properties when it is perturbed 
by a stochastic spatial deformation that can have in turn 
an internal degree of spatial organization, that is, spa- 
tial correlations. Depending on the physical application 
and the context, the perturbation can be seen either as 
a fluctuation due to a physical process or as noise. The 
fundamental question consists in finding how the spatial 
correlation properties of the point process change under 
the effect of the perturbation, and how effective this can 
be in changing the spatial correlations of the point pro- 
cess. 

In this paper we focus our attention on the changes in- 
duced on the two-point spatial correlations of a point 
process by a stochastic displacement field both with 
and without displacement-displacement correlations. We 
work in the hypothesis of statistical independence be- 
tween the point process and the displacement field. The 
exact results presented in this paper can find applica- 
tion in many scientific topics. For instance, in the con- 
text of the so-called n— body cosmological simulations 
, performed to study the problem of "structure forma- 
tion" (e.g., galaxy formation) from the primordial mat- 
ter density field under the effect of the internal gravita- 



tional interaction, point processes are used to represent 
the evolving matter density field. The initial conditions 
of these simulations, representing the primordial density 
field whose spectrum of fluctuations is predicted by the- 
oretical models [l3| , are usually built by applying an ap- 
propriate stochastic displacement field to extremely or- 
dered pre-initial configurations of the point-particle dis- 
tribution [Til Hoi (i.e., either a lattice or a particular 
"glassy" configuration). However in literature the effect 
of the displacement field is described only in approximate 
ways neglecting the contribution of the internal corre- 
lations of the particle system before the application of 
the displacement fie ld J l6j and/or assuming sufficiently 
small displacements |17| using the so-called Zeldovich ap- 
proximation. In the present paper we give the exact de- 
scription of these effects at every spatial scale and for 
any spatially stationary stochastic displacement field. In 
the cosmological context these results can be useful for 
example to understand better the role of the small spa- 
tial scales on the dynamics of the structure formation 

Guam. 

Another important application deals with the prob- 
lem of realizability of point processes with an arbitrary 
a-priori given two-point correlation function satis fying 
the hypothesis of the Wiener-Khinchin theorem [2lll22|. 
This problem is of great importance for instance to study 
the permitted disordered configurations of hard spheres 
systems |24|. While for continuous stationary stochas- 
tic processes (e.g., Gaussian processes) the hypotheses of 
the Wiener-Khinchin theorem give the necessary and suf- 
ficient conditions for the realizability of the process, this 
is not true for point processes. In this case the same hy- 
potheses, adapted to point processes, provide only neces- 
sary conditions for the realizability of the process. There- 
fore finding some limitations in the realizability of point 
processes can be extremely useful for this problem. In 
this context we will show that perturbing a regular lat- 
tice (which can be considered to be the most regular and 
uniform point process) with any stochastic displacement 
field with a continuous spectrum, it is not possible to 
generate a point-particle distribution with an arbitrar- 
ily small degree of disorder with respect to the initial 
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lattice. In fact, we will see that a kind of lower limit 
appears for the degree of disorder injected into an ini- 
tial regular lattice by any displacement field. This lower 
limit is measured by a maximal finite value of the expo- 
nent of the power spectrum of the lattice perturbed by 
the displacement field at small wave numbers. 

Other possible applications come from the study of the 
void distribution and Voronoi tessellation in the super- 
homogeneous class of point processes p^ |. 

The paper and the presentation of the results are or- 
ganized as follows: In Sec. [H] basic statistical properties 
of point processes are briefly presented. This includes 
a classification of all spatially stationary point processes 
in three classes (i.e., essentially Poisson, superhomoge- 
neous, and critical) in terms the asymptotic scaling be- 
havior of the number fluctuations. 

In Sec. IHII we introduce an approximate argument, of- 
ten used in many physical applications, about the effect 
of the displacement field on the two-point correlation 
properties of a point process. In this presentation we 
make clear that this approximation is valid in the limits 
of small displacements and large spatial scales. 

The rigorous treatment of the problem is introduced 
in Sec. ED 



II. BASIC DEFINITIONS 

First of all let us recall some useful definitions 
about stochastic mass density fields. Given a generic 
(discrete or continuous) stochastic mass density field 
p(x) with spatially stationary statistical properties in a 
d— dimensional Euclidean space its average value is de- 
fined by 

<>( x )) = Po , 

where the symbol (...) indicates the ensemble average 
[33|. We will limit our analysis to the case of stochas- 
tic mass fields which can be considered spatially uniform 
(elsewhere homogeneous) at sufficiently large scale. This 
implies that po > 0, excluding in this way the case of frac- 
tal mass distributions for which po — asymptotically 
but with large fluctuations in the conditional density at 
all scales [25| . 

The main correlation properties of the density field are 
given by the covariance function (CF), also called con- 
nected or reduced two-point correlation function: 

C(x - y) - (p(x)p(y)) - pi . 

Another important quantity to characterize the internal 
two-point correlation properties of a stochastic field is the 
so-called power spectrum (PS) 5(k) (also called structure 
factor). It is defined by 

(Mk;£)| 2 ) m 



where 

6 p (k; L) = f f L/2 d d x [p(x) - p ]e- Jkx . (2) 

J J-L/2 

Clearly in the limit L — > +oo Eq. [21 becomes the Fourier 
transform of the density contrast p(x) — po- The quantity 
5(k) measures the net weight of each k— mode to the 
determination of the stochastic process without taking 
into account the phase contribution. 

Due to the spatial stationarity of the stochastic field, 
it is simple to show that S'(k) is equal to the Fourier 
transform of C(x) [2(| : 

S(k) = J d d x C(x)e- 4k - x = FT[C(x)) , (3) 

implying in turn 

c(x) = jhyil ddk s{kyk ' X ^ FT ^ S ^ ■ 

Note that the condition that C(x) vanishes for |x| — > +oo 
implies that k d S(k) -> for |k| -> 0. If the particle 
distribution is also statistically isotropic C(x) depends 
only on x = |x| and S(k) on k — |k|. 

We devote the rest of the paper to the so-called spa- 
tially stationary point processes (SPP), i.e., stochastic 
mass fields consisting of point-particles of unitary mass. 
For this class of systems the microscopic mass density 
p(x) coincides with the microscopic number density n(x) 
which can be written as 

n(x) = ^(x- Xi ), (4) 

i 

where Xj is the spatial position of the i th particle of 
the system, 5(x) is the usual d— dimensional Dirac delta 
function, and the sum is extended to all the particles 
of the system. As aforementioned, we limit the dis- 
cussion to SPP characterized by a well defined average 
number density uq > (i.e., excluding fractal- like par- 
ticle distributions). Due to Eq. 0] and to the fact that 
no > 0, it is simple to find that, for a SPP, the covari- 
ance C(x) = (ti(xq +x)rt(xo)) — 71$ has a diagonal sin- 
gular Dirac delta-like contribution at x = 0. That is, it 
can be written as 

C(x) = n S(x) + l%lh(x) , (5) 

where 7Jq/i(x) is the off-diagonal part measuring the spa- 
tial correlation between number fluctuations in differ- 
ent spatial points, i.e., for x > 0. In general, for truly 
stochastic point processes, this is a rather smooth func- 
tion of x and goes to zero for x — > +oo, but in some 
cases (see below the examples of the "shuffled" lattices) it 
can present also singularities and Dirac delta-like spikes. 
Note that the spatially stationarity and the exclusion of 
a fractal-like behavior imply that there is a finite length 
scale Ao > beyond which fluctuations of the mass con- 
tained in volume of such size become "small" with respect 
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to the average value of the mass itself. Well beyond this 
distance the mass (i.e., number) distribution can be con- 
sidered with good approximation spatially uniform or ho- 
mogeneous and for this reason it is called uniformity or 
homogeneity scale |25j- It is simple to show that in the 
case in which h(x) is sufficiently regular, this scale can 
be defined as the distance Ao such that |/i(x) < f for 
x > \ . From Eqs.|21and|S]we can write the PS of a SPP 
as: 

S{k) = n + nlh(k) , 

where h(k) = FT[h(x)] which in general is a regular 
function decreasing to zero at large k. For point processes 
S(k) in statistics is also called Bartlett's spectrum [27| . 

Since it will be useful to develop the arguments of the 
following sections, we now introduce a brief classification 
fl8| of all the spatially stationary stochastic mass fields 
p(x) (point process or continuous stochastic field) with 
well defined po > in terms of their large scale correla- 
tions and fluctuations: 

(1) If S(k = 0) = c > then C(x) decreases to zero 
at large x faster than x~ d and has a positive integral 
over all space (equal to c), i.e., two-point correlations 
are short range and mainly positive. Moreover, calling 
M(R) = f s , R \ d d x p(x) the mass in a given sphere S{R) 
of radius i?, the average quadratic fluctuation of this 
quantity behaves as (AM 2 (i?)) ~ R d . For this reason 
this class of systems can be called substantially Poisso- 
nian, as on sufficiently large scales the system shows ba- 
sically Poissonian fluctuations. A characteristic physical 
example is given by a homogeneous gas at high temper- 
ature. 

(2) If S(k) ~ k? at small k with -d < [3 < then 
C(x) ~ x~ l3 ~ d at sufficiently large x. That is the sys- 
tem has long range and mainly positive correlations (i.e., 
Jd d xC(x) = +oo). This implies (AM 2 (i?)) - R d ~P 
and for this reason such systems are called superPoisso- 
nian or critical. A physical example in this class is given 
by the density field of a fluid at the critical point of the 
gas-liquid second order phase transition. 

(3) If S(k) ~ fc' 3 with (3 > at small k then we can 
say that at large x the CF C(x) decays faster than x~ d 
and that J d d xC(x) = 0. This means that two-point 
spatial correlations are essentially short range. However 
they are not mainly positive: The condition S(0) = 
indeed implies a precise balance between positive and 
negative two-point correlations. More precisely the rela- 
tion S(0) = can be seen as a condition of geometrical 
order in the spatial organization of the stochastic mass 
fluctuations. The higher is (3, the higher the large scale 
degree of order. As a matter of fact, as shown below, 
in the case of a regular and periodic lattice of particles, 
which is the most ordered particle distribution, one can 
say that S(k) ~ k with (3 -> +oo for k -> 0. For all 
these reasons for this last class of mass fields the name 
superhomogeneous has been proposed |l8l Il9l |25| (else- 
where hyperuniform p^). 



III. AN APPROXIMATE ARGUMENT 

Before entering the detailed and rigorous discussion, 
we give an argument usually used to roughly describe 
the effect of a displacement field on a "sufficiently uni- 
form" mass distribution (for another approximate result 
to this problem going further than the present approx- 
imation see 0|). This argument is based on the fact 
that, if the applied displacements are sufficiently small, 
the mass is conserved "locally". Hence a form of con- 
tinuity equation has to be satisfied. Let us call /9; n (x) 
the initial microscopic density field, and p(x) the same 
quantity after the application of the displacement field 
u(x). By considering the displacements "small" enough, 
we can write the one-step continuity equation: 

p(x) - p in (x) + V • [ j o m (x)u(x)] ~ , (6) 

where u(x) is the displacement performed at point x in 
the given temporal step. The equality is rigorously sat- 
isfied only in the case of infinitesimal displacements. Let 
po > be 

po = (pin(x)) 

It is simple to verify that this average value is not mod- 
ified by the action of the displacement field. If p„(x) 
is "sufficiently uniform" with respect to p(x) , we can ap- 
proximate it in Eq. [S]with a continuous and uniform den- 
sity field pi n (x-) — po, so that Eq. [B|can be rewritten as 

p(x) -p Q + PoV • u(x) ~ . (7) 

By taking the Fourier integral in a cubic volume of size 
L of Eq. [7| and using Eq. |3 we have: 

\5 p (k;L)\ 2 ^pl\k-v(k;L)\ 2 , (8) 

where 

r r L / 2 

v(k;L) = / / Ae- ,lM u(x). 

J J-L/2 

From Eqs. ISJand^ we can say that the PS of the final 
mass distribution is roughly proportional to k 2 times the 
PS of the displacement field (actually the PS of the vector 
field u(x), as its CF, is a d x d matrix and the form of 
such relation can be more complicated). 

This simple result is based on two assumptions: the 
former consists in approximating the initial microscopic 
density with a completely uniform continuous mass field. 
As a consequence we expect that in the exact equations 
describing the effect of the displacement field on a spatial 
distribution of identical particles, there will be terms re- 
lated to the inhomogeneities (e.g., "granularity") of the 
initial particle density. The latter is the fact that, behind 
Eq. there is the assumption of "sufficiently small" dis- 
placements. Then we expect that, when this assumption 
is not valid, the relation Eq. [S] will change drastically. 
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realizations of the displacement field u(x) is defined by 



FIG. 1: The figure presents a pictorial view of the effect of a 
stochastic displacement field to a spatial particle distribution 
in 2d. The particles pass, through the displacements (dashed 
arrows) , from the old positions (black circles) to the new ones 
(empty circles). 



IV. DEFINITION OF THE PROBLEM 

We start by considering a SPP, as defined above, with 
microscopic density rij n (x) given by Eq.0]and with uq > 
0. 

Let us now suppose of introducing a stochastic dis- 
placement field u(x) displacing each particle from its 
initial position. In general this displacement process 
changes the PS of the initial point distribution from 
5i„(k) to a new S(k) (or equivalently the initial CF 
Ci n (x) to a new C(x)). If is the displacement ap- 
plied to the particle i, the position of this particle passes 
from Xi to Xi + Ui (see FigOJ. Therefore the final particle 
density field can be written as 
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N 



n d* 

3 = 1 



/jv(Ul, ...,Uiv)A(ui, ...,Uat) 



(10) 

where fx is the joint probability density function (PDF) 
of the displacements {ui, u^r} applied respectively to 
the set of points {xi, ...,Xjy}- In general, depends 
parametrically on the points positions Xj. In the case of 
a statistically stationary displacement fields, /at depends 
parametrically only on the separation vectors between all 
the couples of the points of the set {xi, ...x^}. A par- 
ticular and very important case is when the set of points 
coincides with the positions occupied by all the particles 
(in which case we call the joint displacement PDF simply 
V({ui} j) of the initial SPP or by the particles of one of 
its subsets. Note, however, that in our hypothesis the 
form of this PDF does not depend on the fact that these 
points are actually occupied by particles. 

Finally, if we have a function of the final (i.e., after 
the introduction of the displacements) microscopic den- 
sity n(x), the ensemble average over all the possible final 
configurations of the particle distribution is given by av- 
eraging over all the possible displacements as in Eq. 1101 
fixing the initial particle density ni„(x), and then over 
all the possible initial particle configurations (...). This 
is due to the fact that the ensemble of the all possible fi- 
nal particle configurations is found by considering all the 
possible initial configurations, and for each of these all 
the final configurations obtained by applying the ensem- 
ble of the displacement fields. However, if, as we suppose 
here, the displacement field is statistically independent 
on the initial particle distribution, the order of these two 
averages is arbitrary. For instance, in this case, the CF 
of the "displaced" particle distribution is expressed by 



C(x) = ( n(x + x)n(x ) 



n-(x , 



(11) 



n(x) 



6(x 



(9) 



with an arbitrary order of the two averages (...) and 



A stochastic displacement field u(x) can be seen as a 
continuous stochastic vector field. We assume that this 
field is spatially stationary in the statistical sense, i.e., it 
is characterized by the invariance of the statistical prop- 
erties for any spatial translation. We can think to "at- 
tach" a displacement vector u(x) to each spatial point 
x, even though it acts on the mass density only if x is 
occupied by a particle. In what follows, we assume that 
the displacement field is statistically independent of the 
realization rii n (x) of the initial particle density, i.e., the 
probability of having a given realization u(x) of the dis- 
placement field is independent on the realization of the 
initial particle distribution. 

Let us consider a function A only of the displacements 
{ui, ujv} applied respectively to a set of spatial points 
{xi, xjv}. The average of this quantity over all the 



V. EXACT RESULTS FOR ONE AND 
TWO-POINT STATISTICAL PROPERTIES OF 
THE PARTICLE DISTRIBUTION 

The aim of this section is to relate the one and two- 
point correlation properties of the "displaced" particle 
distribution to its initial ones and to those of the applied 
displacement field by finding exact relations going beyond 
the approximation given by Eq. 1111 As aforementioned, 
the discussion will be limited to the case of spatially sta- 
tionary stochastic displacement fields and initial particle 
distributions. In this way also the final particle distribu- 
tion will be spatially stationary. 

We will start by evaluating, through Eq. the aver- 
age mass density (n(x)V The next step will consist in 
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finding the transformation equation for the PS S(k) (or 
equivalently the CF C(x)). 

Since the displacement process does not create or de- 
stroy any particle and is statistically stationary, the av- 
erage mass density stays equal to the initial one no: 

n(x)^> = n ■ 

This can be also proved by direct calculation using 
(nj n (x)) = no. First of all we note that Eq.[5]is a sum of 
single displacement terms. Therefore, in order to evalu- 
ate the displacement average n(x), we need only to know 
the one displacement PDF /i(u) and not the complete 
joint PDF "P({u,;}) for all the system particles. In our hy- 
pothesis of spatial stationarity /i (u) does not depend on 
the point of application of the displacement, and we recall 
it p(u) for simplicity (p(u) is obtained from V({ui}) by 
integrating out all the displacements with the exception 
of one). We can then write 

n(x) = / d d u z p(ui)S{x -Xi- Uj) = J^p(x - Xj) . 

i i 

By taking the average (...) over the initial particle con- 
figurations, we finally have 



= y <2<W)£>(y-x + x; 
n / d d yp(y) = n , 



where we have used the statistical spatial stationarity of 
n* n (x) (i.e., (n in (x)) = (n in (y - x)) = no) and the nor- 
malization condition of the one-displacement PDF p(u). 

We can now face the problem of calculating the new 
CF C(x) and the new PS 5(k). The key point is to 



evaluate the average (n(x)n(y) ). Since the product 



n(x)n(y) 



i,3 



6(x - Xj - Ui)5(y - Xj - Uj 



is a sum of terms containing either one (for i = j) or 
two (for i 7^ j) different displacements, we do not need 
to know the complete joint PDF ^({uj}), but only the 
joint two-displacement PDF /2(u, v) which is obtained 
from V({ui}) by integrating out all the displacements 
but two. In general /2(u, v) will depend parametrically 
on the coordinates of the two points of application of the 
displacements. Assuming the hypothesis of a spatially 
stationary displacement field, /2(u, v) depends paramet- 
rically only on the separation vector x between these two 
points. For this reason we recall it / 2 (u, v) = /(u, v;x) 
putting in explicit evidence this dependence. Note that 
the function /(u, v;x) carries much more information 
than the simple knowledge of the average displacement 
u = U, and the two-displacement correlation matrix of 
elements 



with fi, v = 1, ...,d, where is the [i th component of 
the displacement vector u. In fact (x — y) is only the 
average value of («W(x) - U^) (u M (y) - C/ M ) calcu- 
lated with the PDF /(u, v; x) itself, while the knowledge 
of /(u,v;x) include all information about all the higher 
moments of two-displacements. 

The joint two-displacement PDF /(u, v; x) satisfies the 
following limit conditions in x: 



/(u,v;0) = <5(u-v)p(u) 
lim /(u,v;x) = p(u)p(v) 



(13) 
(14) 



The former equation is trivial, while the latter estab- 
lishes simply that the correlation between two different 
displacements must go to zero if the distance between the 
two points of application goes to infinity. 

First of all let us evaluate the average of n(x)n(y) over 
the displacements. By direct integration one obtains 



"(x)n(y) = (15) 
^2 d d u i d d u J [f(u i ,u :j ;x l3 )S(-x - Xj - Uj) x 



h3 



S(y - Xj - u0] = ^ /(x - x*,y - Xjjxy) 



where Xy- = x$ — Xj. Note that the first limit condition in 
Ea. ll3l pcrmits to perform the average without separating 
the diagonal contribution i = j from the off-diagonal part 
ij^jof the double sum in Ea. ll5l bv averaging separately 
the former using the one-displacement PDF f>(uj) and the 
latter through the two-displacement PDF /(uj, u^-; x,j) 
with i ^ j. 

The next step is to evaluate the average (...) on the en- 
semble of initial particle configurations. For this purpose 
note that the ensemble average of a function of the initial 
configuration of the form J2 i ■ ijjfe, Xj), where y) is 
a generic two-point function, can be written as 



5>( Xl ,x,)^ (16) 
d d xd d y V(x,y) - *i)*(y - = 

d d xd d yi;(^y)(j2S(x-^)6(y-^ = 
d d xd d y (n in (-x)n in (y)) ^(x,y) , 



G^(x - y) = («(*») ( x ) - UM) (uM(y) - UM) (12) 



where, by definition, we have 

(n in (-x.)n in (y)) = n„ + C m (x - y) . (17) 

Note that the diagonal part n 5(x) of the connected two- 
point correlation function Ci n (x) takes correctly into ac- 
count the diagonal term i — j of the sum of Eq. 1161 
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By applying Eas.lTBlandlTTIto Eq.^J it is possible to 
write: 

(n(x)n(y)) = j J d d x'd d y' [n 2 + C m (x' - y')] 
x/(x-x',y-y';x'-y') (18) 

It can be convenient to rewrite Eq. 1181 by separating the 
two terms coming respectively from the diagonal and the 
off-diagonal parts of Cj„(x), that is by writing, as in 
Eq.El Ci n (x) = n 5(-x) + nlh in {yL): 



|2lj')- We can now apply Eq.[2]to Eq.^|in order to find 
the two-point correlation function of the final system: 



n ( x My)) = n S(x - y) 



(19) 



.„ , , d d x'd d y'll + h in (x'-y')] 
x/(x-x',y-y';x'-y'). 

We are now able to write the new CF C(x) of the final 
particle distribution which is defined, as usual, by 



C(x) = (n(x + x)n(xo)^ - n 2 , . 



(20) 



Note that from Eq.^lthe diagonal part of C(x) remains 
equal to that of Ci„(x) (i.e., no<5(x)) as expected. 

In order to write the transformation equation for the 
PS, which is the most useful in many applications [T3 |. 
we start from the simplest case of uncorrelated displace- 
ments; then we will come back to the general case for 
general considerations and some paradigmatic examples. 



VI. INDEPENDENT DISPLACEMENTS 

We now consider the case in which the displacement 
applied to a given spatial point is statistically indepen- 
dent of the displacement applied to any other point. 
Therefore the statistics of the stochastic displacement 
field is completely determined by the knowledge of the 
reduced one-displacement PDF p(u), and the joint PDF 
of n displacements Ui,U2, ...,u„ in n different points of 
the space factorizes as follows: 



/„(ui,u 2 , ...,11,, 



In particular for the two-displacement PDF, we can write 

{5(u — v)p(u) for x = 
(21) 
p(u)p(v) for x 7^ 

Note that the lack of displacement-displacement correla- 
tions implies a discontinuity of /(u,v;x) at x = 0. As 
shown below this does not happen for truly continuous 
correlated stochastic displacement fields (i.e., belonging 
to the class of continuous stationary stochastic processes 



n(x)n(y)^ = n 2 , + n 5(x - y) 



(22) 



d d x'd d y' p(x - x')/i m (x' - y>(y - y') 



Since in this case /(u, v;x) is discontinuous at x = 0, it 
has been important to separate the contributions of the 
diagonal and the off-diagonal part of Cm (x) in Eq. 
In fact, in the hypothesis of uncorrelated displacements, 
any element of the connected two-displacement correla- 
tion matrix G A1!/ (x), given in Eq. ^1 vanishes for any 
x ^ while G M1/ (0) = S^g 2 where g 2 > is the single 
displacement variance. That is, G Atl/ (x) is discontinuous 
at x = 0. As aforementioned, this is a very particular 
case, as in truly correlated continuous stationary stochas- 
tic processes it is continuous everywhere pl| . 

At this point, by remembering that the PS 5(k) = 
FT[C(x)], with C(x) given by Eq. 1201 we can Fourier 
transform Eq.[221to obtain the following local relation in 
k space: 

S(k) - no(l - |p(k)| 2 ) + |p(k)| 2 5 m (k) , (23) 

where p(k) is the characteristic function of the one- 
displacement PDF 

p(k) = FT[p(u)] , 

and where we have used no-FT[/i iT j(x)] = Sj n (k) — no- 
Note that by definition p(0) = 1. 

Equation 1231 gives the relation between the PSs of the 
point-particle configurations before and after the appli- 
cation of the uncorrelated displacements field. First of 
all let us analyze the notable case in which the initial 
point-particle distribution is the statistically stationary 
and isotropic Poisson one, i.e., that case in which there 
is no correlation between the init ial p ositions of different 
particles. It is simple to show [l8ll2fJ | that the initial den- 
sity CF is simply Cj n (x) = no<5(x) (i.e., it has only the 
diagonal part). This means that Si„(k) = no which, in 
view of Eq. |23 implies 5(k) = n too, regardless to the 
form of p(u). That is, the particle distribution remains 
Poissonian after the application of any random and un- 
correlated displacement field. This can be considered as 
a formulation of the so-called theorem of Bartlett [3lJ (for 
the behavior of a Poisson point process under a determin- 
istic displacement field see |32Tp. This property is quite 
easy to understand: the displacement field has no spatial 
correlation, hence it tends to randomize the particle dis- 
tribution, but the Poisson particle distribution is already 
the "most random possible" SPP. This is further clarified 
by observing that uncorrelated displacements cannot in- 
crease the spatial correlations in the particle distribution. 
We have just shown that the stationary Poisson point 
process of average density no > is the "fixed point" of 
the transformation given by Eq. 1231 We show now that 
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this fixed point (i.e., the stationary Poisson point pro- 
cess) is also attractive. That is, we show that, starting 
from an arbitrary stationary point process with density 
no and initial PS Si n (k), by applying successive stochas- 
tic uncorrelated displacements to all the particles taken 
from the same PDF p(u) with no correlation at each step, 
the PS flows toward Soo(k) = no, i.e., the Poisson one. 
It is simple to show that after m consecutive applica- 
tions of such displacement field the PS 5 m (k) satisfies 
the relation: 

S m (k) = no(l - |p(k)| 2m ) + \p(k)\ 2m S m (k) . 

The previous equation is simply obtained by the fact that 
the characteristic function of the sum of m independent 
random vectors extracted from the same PDF p(u) is 
simply [jpfk)!" 1 . Since p(u) is a PDF, a part from partic- 
ular cases in general \p(k)\ < I for any k > 0. This 
implies that 5 m (k) converges exponentially fast to nn for 
each k, i.e., to the PS of the Poisson point process |3q . 

One has also to notice that the right hand side of Eg. 1231 
is the sum of two terms: the former is proportional to the 
average density of particles no, that we can call granular- 
ity term, and independent of the initial PS of the particle 
distribution, while the latter depends on no, only through 
Sin(k) which satisfies the condition Si n (k — > oo) = no 
because of the diagonal term of Cj n (x) and of the fact 
that the off-diagonal part must be integrable at small x. 

If the initial point process is statistically isotropic as 
well as stationary, then Cj„(x) is a function only of x = 
|x| and Sj n (k) only of k = |k|. Furthermore, if also the 
displacement field is statistically isotropic, then = 
for each fj, = 1, d, and p(u) depends only on u = |u|. 
This implies that also S(k) depends only on k (and C(x) 
on x). 

A. Small k expansion and large scale behavior - I 

In many applications, such as cosmological studies 0, 
it is particularly important to analyze the behavior of 
S(k) at small k, that is at large spatial scales. For this 
reason we now study the asymptotic behavior of Eq. 1231 
for k — > 0. We limit the discussion to the case in which 
both the point process generating the initial particle dis- 
tribution and the stochastic displacement field are statis- 
tically isotropic. As seen above, this hypothesis implies 
S ln (k) = S in (k), p(u) = p(u), and S(k) = S(k). 

The first step consists in studying the small k behavior 
of the characteristic function p(k). By definition we have 




and then p(0) = 1. As p(u) — p(u), then p(k) = p(k). 
Let us suppose that at sufficiently large u we have p(u) = 
Au~ a + o(u~ a ), where a > d as p(u) must be by defi- 
nition integrable over all the space (a — > +oo includes 
exponential-like or faster decay at large u). Using this 



property and the definition of p(k), it is simple to show 
that to the lowest order in k larger than zero, one has 

( /3 = ot-d if0<a-d<2 
p(k) = 1 - Bk 13 with I 

{ (3 = 2 ifa > d + 2 

(24) 

where B > 0. Moreover if a > d + 2 then B = ||. 
Instead, if d < a < d + 2, then u 2 diverges and 

B = aJ d d xx~ a (1 - e - lxcos0 ) , (25) 

where 9 is the angle between x and any one of the co- 
ordinate axes. Note that in any case < (3 < 2. This 
implies that for k <C B 1 !^ 

1 - |p(fc)| 2 ~ IB kP . 

On the other hand, as seen in Sec.^ in order to have the 
initial SPP well defined, it is necessary that k d Si n (k) — > 
for k — > 0, that is practically St n (k) ~ k 1 at small k with 
7 > -d. 

We can draw the following conclusions for the small k 
behavior of S(k) in Ea.l23l 

1. If — d < 7 < (with as seen above < (3 < 2), 
S(k) ~ k 1 similarly to Si n {k), and its small k am- 
plitude is independent of the displacement field. 
This means that, in this case, finite uncorrelated 
displacements cannot destroy the persistence of 
correlations already present in the system. In par- 
ticular it is important to note that if 7 < (i.e., 
long range density-density correlations in the initial 
particle configurations) no uncorrelated displace- 
ment field is able to affect the large scale correlation 
properties of the initial system. 

2. On the contrary, if < (3 < 7, then the small k 
behavior is completely determined by the displace- 
ment field, resulting in S(k) ~ 2Bk^. As shown in 
Sec. m a point-particle distribution having 7 > 
is called superhomogeneous showing a sort of long 
range order, characterized by subPoissonian mass 
fluctuations at large scales, which in the present 
case is partially destroyed by the Poissonian noise 
introduced by the displacement field. In this re- 
spect, note that if u 2 = +00 then (3 < 2. Conse- 
quently, it is much simpler to obtain the condition 
7 > j3, i.e., the randomization of the system in- 
troduced by the uncorrelated displacement field is 
more effective. 

3. If (3 = 7 both the long wave length modes of the ini- 
tial particle configuration and of the displacement 
field determine of the small k power spectrum of the 
final system. In particular the exponent is equal to 
the initial one while its amplitude increases. This 
indicates that the initial and the final systems have 



FIG. 2: The projection on the xy plane of a 3d shuffled lattice 
of 16 3 particles in the unitary volume is represented. In the 
present case each particle is randomly displaced from its ini- 
tial lattice position inside a finite cubic box centered around 
the lattice point and of side equal to one fifth of the lattice 
spacing. The displacement applied to each particle is statis- 
tically independent of the displacements applied to the other 
particles. 

the same kind of long range order of density fluc- 
tuations and the same mass-length scaling relation 
for the large scale mass fluctuations, but with an 
increase of their amplitude in the second case. 

B. The shuffled lattice with uncorrelated 
displacements 

In this subsection we present a simple, but impor- 
tant example of application of uncorrelated stochastic 
displacement fields: the random shuffling of a regular 
lattice of particles. Its importance relies on the fact that 
a perturbed lattice is often used as the initial condition 
for many dynamical applications as, for instance, the al- 
ready mentioned cosmological n— body simulations |l4j . 
and bio- metrical studies In this section we present 
the simplest example of a lattice with stochastic displace- 
ment perturbations. In Fig. [21 the projection on the xy 
plane of a three dimensional shuffled lattice is given. 

It is simple to show that for a distribution of point 
particles of unitary mass occupying the sites of a regular 
cubic lattice, the PS is [l| 

S in (k) = (2n) d n 2 J2 S (k - H) , (26) 

H#0 

where the sum runs over all the sites H of the recipro- 
cal lattice P, with the exception of the origin 0. We 



FIG. 3: Power spectrum of a Id shuffled lattice with un- 
correlated displacements with finite variance obtained both 
from numerical simulations and from the theoretical predic- 
tion given by Eg. 1271 and showing perfect agreement. In par- 
ticular the two PSs refer to a regular chain of particles and 
unitary lattice spacing perturbed by the displacement field, 
through which each particle is randomly displaced in a point 
of the segment centered around its initial lattice position and 
of length a = 1/50 (i.e., p(u) = 9{a/2-u)6(u+a/2)/a). Note 
that at small k the PS S(k) scales as k 2 as the displacement 
variance is finite. Moreover the Bragg peaks (whose ampli- 
tude here has been normalized for pictorial reasons) are per- 
fectly modulated by weights proportional to |p(fc)| 2 (which in 
this case is given by p(k) = 2 sin 2 (ka/ 2) /ka). Finally the PS 
at large k converges correctly to the average number density 
n = 1. 

remind that, if the particles occupy the sites of a cubic 
lattice with lattice spacing I, then each component of H 
is a positive or negative integer multiple of 2n/l and the 
average density of points is uq = l~ d . 

We can now apply Eq. [23 in order to find the final PS 
S(k) after the random shuffling (i.e., the application of 
the random displacement field): 

5(k)-n (l-|p(k)| 2 )+(27r) d n^ |p(H)| 2 <5 (k - H) , 

H^O 

(27) 

We now stress two important aspects of Eq. ITH 

1. The random shuffling in general does not erase com- 
pletely the presence of the so-called Bragg peaks 
(i.e., the sum of delta functions), but only modu- 
lates their amplitude and adds a continuous con- 
tribution typical of fully stochastic point distri- 
butions. The complete cancellation of the Bragg 
peaks contribution to S(k) is possible only in the 
very particular case in which p(H) = for every 
reciprocal lattice vector with the exception of 0. 



2. around k = (more precisely in the so-called 
first Brillouin zone |2j of the reciprocal lattice) 
<Si n (k) = at any order of the Taylor expansion 
(we can loosely say that Si n (\c) ~ fc + °° around 
k = 0). Consequently, as clear from Eq. [23 in 
this region 5(k) is determined by only the behav- 
ior of the displacement characteristic function p(k) . 
As shown above, if the displacement field is sta- 
tistically isotropic, p(k) = p(k). Therefore, even 
though the lattice is strictly anisotropic, the shuf- 
fled one has isotropic mass fluctuations at large 
scales. This implies that on sufficiently large scales 
the scaling exponent of the fluctuations of the num- 
ber of particles contained in a given volume of lin- 
ear size R do not depend strongly on the shape 
of the volume itself. This is not true, instead, for 
a deterministic cubic lattice for which the scaling 
behavior of these fluctuations changes in passing 
from a spherical to a cubic volume with the same 
symmetry of the lattice. 

Since in the first Brillouin zone S(k) is completely de- 
termined by p(k), the asymptotic behavior at small k 
of S(k) can be derived by Eq. [21 In particular if the 
variance u 2 of the displacement field is finite, we find 
S(k) ~ k 2 for k — » 0, independently on the particular 
form of p(u) (see Fig. [31 • This is a case of universal be- 
havior for all the PDF p(u) with sufficiently fast decay 
at large u. 

Instead, in the case in which u 2 diverges, i.e., d < a < 
d + 2, this universality is lost, having 5(k) ~ k" with 
f3 = a — d 7 with a one-to-one correspondence between 
the exponents a and as shown in the previous section 
(see Fig. @J. A similar case of universality is found in 
random walks with independent steps [28J. In fact if the 
variance of the steps is finite (ordinary random walks) 
the average quadratic distance (/S.x 2 (t)) reached by the 
walker after a large number t of steps satisfies the scal- 
ing relation (Air 2 (£)) ~ t independently of the precise 
functional form of PDF of the single step. On the other 
hand, if the single step variance goes to infinity (Levy 
flights) this is no more true, (Ax 2 (t)) being infinite, and 
the PDF of Ax(t) at sufficiently large t having a power 
law tail with an exponent in a one-to-one correspondence 
with that characterizing the tail of the PDF of a single 
step. A similar transition from a universal scaling behav- 
ior of the fluctuations to a non-universal one has been 
found also in more complex fragmentation problems [29j| . 



VII. CORRELATED DISPLACEMENTS 

Let us now go back to Eq. ^3 and consider the general 
case of a stationary stochastic displacement field with 
spatial correlations. In this case /(u, v;x) cannot be 
factorized as in Eq. [21 for x > 0. 

In order to write the equation of transformation of the 
PS let us remind the basic relation between the CF and 
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FIG. 4: Shuffled lattice of the same kind of that in Fig. [|1 
but now with p(u) = {a/n)/{u 2 + a 2 ) where a = 1/50, i.e., 
with unlimited displacements and logarithmically diverging 
u. Again the agreement between the numerical simulation 
results and the theoretical prediction Eq. l27l is excellent. Note 
in particular that S(k) ~ k at small k and that the remaining 
Bragg peaks, superimposed to the continuous contribution to 
the PS, are well modulated by weights proportional to \p(k)\ 2 
with p(k) = exp(— a\k\). Finally the PS at large k converges 
to average number density no = 1. 

the PS of a spatially stationary stochastic process: 

J J d d xd d y e- l(k x+q y) C(x - y) = (2ir) d 6(k + q)S(k) . 

(28) 

We also recall that by definition: 

(n(x)n(y) \ =n 2 a + C(x - y) . 

Furthermore, we define the functions /(ki,k2;x) and 
F(ki, k2j q) respectively by the following FTs: 

/(k 1) k 2 ;x)= J J d d ud d v e^ klU+k2V )/(u,v;x), 

(29) 

F(ki,k 2 ; q) = J d d x e-^/Cki, k 2 ; x) . (30) 

The function / , (k 1 ,k 2 ;x) is simply the characteristic 
function of the joint two-displacement PDF. By defini- 
tion /(k!,k 2 ;x) satisfies the following limit conditions 

/(0, 0;x) = 1 for any x, 

and 

/(0, k; x) = /(k, 0; x) = p(k) for any x > . 
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By using Eg. 1191 and Eas. 1281501 we can write: 
S(k) = n (l - j^yl J d d qF(k, -k; q)) + (31) 
d<W ikx /(k, -k; x) (n 2 + C m (x)) - (2n) d n 2 6(k) . 



Note that J d d qF(k, — k; q) must be carefully han- 
dled. In fact, due to the properties of the inversion of 
the Fourier transform, it cannot be substituted directly 
by /(k, — k; 0) if /(u, v; x) is discontinuous at x = and 
continuous anywhere else as in the case of uncorrelated 
displacements presented above. Thus it must be under- 
stood as 



1 



(27r) d 



d d q F(k, -k; q) = lim /(k, -k; x) 



(32) 



More specifically, if /(u, v; x) is continuous at x = 0, i.e., 
lim /(u,v;x) =/(u,v;0) = S(u - v)p(u) , 

x— >0 



then wc have 



y Af(k,-k;q) = l. 



(2 

This condition is valid in all the cases in which the 
stochastic displacement field is a real continuous corre- 
lated stochastic process (see below the Gaussian exam- 
ple). In fact, in this case it is possible to prove a theorem 
[21| stating that the two-displacement correlation func- 
tion is continuous everywhere, being equal to the one- 
displacement variance u 2 — u 2 at x = 0. Instead, in the 
case of an uncorrelated stochastic displacement field this 
is no more true (it is not a continuous stochastic pro- 
cess), and, /(u, v;x) is discontinuous at x = as shown 
by Eq. This together with Eq. [321 gives for this case 



1 



(27r) d 



d d gF(k,-k;q)Hp(k)| 



With this prescriptions it is simple to recover Ea. l23l from 
Eq. El in the case of uncorrelated displacements. 

Instead, in the present case of a stationary correlated 
continuous stochastic displacement field, Eq. 1311 can be 
rewritten as 



S(k) = 



7>,-k;x) [n 2 + C in (x)] 



(33) 



- (2n) d n 2 S(k), 
or equivalently 

5(k) = n§F(k,-k;k) + (34) 
J d d qF(k,~k;q)S m (k- q )-(2TT) d n 2 S(k). 

Equations 1331 can be further simplified by noticing that 
in the case of spatial statistical stationarity of both the 



initial point-particle distribution and of the displacement 
field, the PS of the final particle distribution will not 
depend separately on the couple of displacements u and v 
applied at two points separated by the separation vector 
x, but on the relative displacement w = u v. In fact 
let us call 0(w; x) the PDF that two points, separated by 
the separation vector x, perform a relative displacement 
w. By definition 

</>(w;x) = y | AA/(u,v;x)i(w-u + v). (35) 

If we take the FT of Eq. 123 with respect to w, we have 

0(k;x) = /(k,-k;x), 

where <f>(k; x) = FT[(f>(w; x)]. Therefore Eq. 1231 can be 
rewritten again 

S(k) = J d d xe- lkx j>(k;x) [n 2 + C ln {*)\~(2iT) d n 2 5{k) . 

(36) 

By using the PDF </>(w;x), Ea. 1181 can be rewritten in 
a very intuitive form. This is done by noticing that 
for a generic stationary point process n(x) the quantity 
(n(x)n(y)) /hq (where (...) is as usual the average over 
the considered ensemble of configurations n(x)) repre- 
sents the average conditional density 25] of particles seen 
by a generic particle of the system at a separation x — y 
from it. By calling r in (x - y) = (n in (x)n in (y)) /n and 

ry(x — y) = (n(x)n(y)\ /no the conditional densities 



respectively before and after the application of the dis- 
placement field, Eq. ^1 transforms into: 

r / (x) = Jd d x' r m (x')0(x-x';x'). 

Very often the CF is defined in a dimensionless way 
dividing C(x) by rip, i.e., it is redefined as 



C(x) _ «5(x) 



"o 



h(x) 



which we call normalized covariance function (NCF). 
Consequently, also the PS is redefined as the FT of £(x): 



P(k) = FT[£(x)] = 



S(k) 



1 

n 



h(k). 



It is simple to verify that if all the particles of the dis- 
tribution have the same non-unitary mass m > 0, i.e., 
if the microscopic mass density is p(x) = m n(x) = 
TO Si^( x — x i) w ith po = mno, the NCF £(x) of the 
microscopic mass density p(x) is equal to that of the mi- 
croscopic number density n(x) and then is independent of 
to. For these rescaled quantities Eq. 1361 can be rewritten 
as 

P(k) = J d rf x e - ik x ^(k;x) [l + &„(x)] - (27r) d <5(k) , 

(37) 
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It is important to note that, while £j n ( x ) depends on 
rip at least through its diagonal part 5(x)/no, the two- 
displacement /(u, v; x), and therefore </>(k; x), in our hy- 
pothesis, are in general supposed not to be (unless for 
particular choices of its momenta). Therefore, differently 
from the case of uncorrelated displacements, both Ecis.1551 
and 1351 can be divided into two parts, one dependent on 
no, and the other independent of it. This is a very im- 
portant point because spatial distributions of particles 
with equal masses m are often used in numerical simula- 
tions as a discrete representation of continuous stochastic 
mass density fields. Consequently, in Eq. OS and Ea. 1351 
there is a part depending on the discretization process 
and another part independent of it. This aspect is par- 
ticularly important in the context of the gravitational 
n— body simulations in which, as above mentioned, the 
matter density field is usually represented, for numerical 
reasons, through a more or less dense distribution of par- 
ticles with the same mass m jl4j so that mriQ = po with 
po put equal to the average mass density of the continu- 
ous model. 



A. Small k expansion and large scale behavior - II 

Starting from Eq. 1371 one can write a simple formula 
for the small k (i.e., large scale) behavior of the final PS 
of the particle distribution after the application of the 
displacement field. First of all let us study the small k 
behavior of </>(k;x) = f(k, — k;x). Since it is the FT 
of </>(w;x), if G^ U (Q) < +oo (implying that the average 
value of w 2 is finite at any x) , we can write 



0(k; x) = 1 — ik. ■ w(x) — 



[k • w(? 



+ o(fc 2 ), (38) 



where w(x) is the relative displacement between two par- 
ticles ini tially separated by a vector distance x. Let us 
suppose w(x) = which is automatic in the case of 
statistical invariance for space inversion or any rotation 
(isotropy). Moreover, using the definition of G A1J y(x) we 
can write 



[k • w(x)] 2 = 2^ fcWfcW [G M „(0) - G^(x)] 



where we have supposed also G M „(x) = G M „(— x). Using 
this expression and Eq. 1381 in Eq. 1371 we obtain at small 
enough k 



P(k)~P in (k) 

l.d 



l-^feMftMG^(O) 
d d q 



G M „(k) 



(27T) 



(39) 

d G^(q)Pi n (k-q) , 



where G M ^(k) = FT [G M1/ (x)] is the power spectrum ma- 
trix of the displacement field. Note that, since G M „(0) = 



/ (2ir) d ^WW) Eq. OH can be also rewritten as: 



P(k)~P in (k)- 

(27T)* " 



l.d 

E 



kMkM {cw(k) 



(q) [P m (k-q)-P m (k)] 



which is useful in particular in the case in which the 
initial particle configuration is the stationary Poisson one 
for which Pj n (k) = 1/uq and the last term vanishes. 

Depending on the small k properties of P; n (k) and 
G M i,(k), the small k behavior of P(k) will change. Note 
that Eq.[2Hlis valid only if G^„(0) < +oo. In the case in 
which G M „(0) = +oo then the singular part of the small 
k expansion of 0(k; x) has to be considered in a similar 
way to the case of uncorrelated displacements given by 
E as. 1241 and!25l A very particular and important case of 
Eq.[3§]is given when perpendicular displacements arc not 
correlated at any x. This means that G Ai;y (x) = <5 /Jl/ G(x) 
and G M „(k) = <5 M „G(k) (with 5^ the Kronecher delta). 
If the displacement field is also isotropic, G(k) depends 
only on k (and consequently G(x) on x), so that the 
PS matrix with elements G Mt ,(k) is invariant under any 
spatial rotation. 

In this case Eq. 1391 can be rewritten as 



P(k) ~ P,, 



G(k) 



(k) [1 - G(0)fc 2 
d d q 



(40) 



(2tt) c 



G(q)P m (k-q) 



Similarly to Eq. also Eq. 0HI can be re-expressed as 
follows: 



P(k)~P in (k) 
d d q 



k 2 |G(k) 
d G(q) [P m (k-q)-P m (k)] 



(41) 



Equations and El are very useful because to show 
very clearly how effective a displacement field can be in 
"injecting" large scale correlations into a given particle 
distribution. As better clarified below, a central role is 
played by the spatial dimension d. Let us suppose that 
at small k we have Pm(k) ~ k a and that G(k) ~ k^ (as 
already shown a, f3 > —d). Note that 

d d q 



( 7 G(q) P m (k - q) = J d d x e- 4k - x G(x)£ m (x) , 

and £m(x) goes to zero at large x. Hence, when /3 < 
0, in Eq. 0U| the term fc 2 G(k) is more effective than 

^ 2 / (27r) d ^(cO-PinXk ~ <l) (i- e -j the latter has an expo- 
nent larger than j3 + 2). For j3 — 0, apart from the 
very particular case in which J ^ d Jj d G(q)Pi„(q) = 0, 
both terms are of the same order at small k, i.e., pro- 
portional to k 2 . Instead for f3 > 0, apart again the pre- 
vious particular choices of G(k) in relation to Pi„(fc), 
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fc 2 /(^G(q)F m (k - q) ~ fc 2 prevails on fc 2 G(k) ~ 

k@ +2 at sufficiently small fc. 
We can therefore conclude that 

• /3<0:ifa</3 + 2 the displacement field is ineffec- 
tive in changing large scale two-point correlations 
between particles. In fact the small fc leading term 
of the final PS -P(k) is the initial one P,- n (k). For 
what concerns perturbations to this leading term, 
it is simple to show that, if a < (3 the main pertur- 
bation to the leading term due to the displacements 
is -G(0)fc 2 P m (k) ~ k a+2 , while if/3<a<^ + 2 
it is fc 2 G(k) ~ k@ +2 (if a = (3 the two terms are of 
the same order). 

If instead a > (3 + 2, the displacement field com- 
pletely determines the new large scale correla- 
tion properties of the particle distribution being 
fc 2 G(k) ~ fc /3+2 now the leading term of P(k). 
Since (3 > —d, the limit "most critical" behavior 
of P(k) which can be reached is k 2 ~ d . Note that 
for d > 3 long range non-integrable and mainly 
positive (i.e., critical) two-point correlations can be 
introduced in the system by the action only of dis- 
placements having finite variance G(0). Accord- 
ingly the (unreachable) limit decaying behavior for 
the NCF at large x is given by £(x) ~ x~ 2 which 
is not integrable and long range. For example in 
d = 3 one can start from a completely uncorrelated 
Poisson particle distribution, and, applying finite 
displacements (i.e., with a finite average squared 
value) but with long range correlations, to obtain 
a particle distribution with a covariance function 
similar to those found at the critical point of a sec- 
ond order phase transition. For d < 2 the exponent 
2 — d > 0, hence £(x) decays faster than x~ d and 
correlations are not "critical", but integrable (i.e., 
for our purpose, short range). 

Finally, for — 2 < j3 < and a > + 2 and any d, 
even though the leading term of -P(k) is due to the 
displacement field which is characterized by long 
range correlations, the final particle distribution re- 
mains superhomogeneous (see Sec. ITTf) . 

• (3 = 0. If a < 2 the leading term of the 
final PS P(k) is again P in (k). The displace- 
ment field introduces only higher order pertur- 
bations: (1) -G(0)P m (k)fc 2 - k a+2 , and (2) 



G(k) + J 



d a q 



fc 2 . The 



ptG(q)Pi„(k-( 

former is the most important perturbative term for 
a < ( "critical" initial condition) , while the latter 
is when < a < 2. For a = (1) and (2) are in 
general of the same order and both contribute to 
the main perturbation to Pj n (k). However, if the 
initial particle configuration is a stationary Pois- 
son one, Pin(k) = 1/no for any k. Consequently, 
as shown explicitly by Eq. 1411 the main perturba- 
tion to Pi n (k) is only fc 2 G(k), as the other terms in 



Eg. 1411 cancel one each other. We show this better 
in the following through the example of the Gaus- 
sian displacement field. 

Instead, if a > 2 again the small fc behavior of 
the final PS is completely different from the ini- 
tial one and is determined by the displacement 
field. In fact now the leading term becomes 



fc 



G'(k) + / I ^G(q)P in (k-q)j ~ , 

ever the system after the action of the displace- 
ments is still in the superhomogeneous class of 
point processes in any d. 



How- 



(3 > 0. If / t|4j G(q)P in (q) > 0, this case is simi- 
lar to the case (3 = with the difference that when 
a > 2 the main term introduced by the displace- 
ment field is only fc 2 / ^ G ( q ) p m ( k _ q ) ^ fc 2 . 

However if J , 2 h G(q) Pi»(q) = 0, which corre- 
sponds to a peculiar choice of the dependence of 
G(q) on Pj n (q), the term fc 2 G(k) can again be im- 
portant. This case will be analyzed in more detail 
in the next section when the case of a shuffled lat- 
tice with correlated Gaussian displacements is an- 
alyzed. For the moment we notice only that this 
case is quite difficult to be realized because of the 
condition given by the Wiener-Khinchin theorem 
|22^ which states that both Pj n (k) and G(k) must 
be nonnegative at any k. This means that it can 
be realized only if P^ n (k) is zero where G(k) is not 
and vice-versa, i.e., non-overlapping supports. 

Finally, note that when the initial particle 
configuration is a regular lattice, the term 
P„(k) [l - G(0)fc 2 ] in Eq. HJ is identically zero 
around fc = 0. Consequently, the small fc PS of the 
final configuration is determined by the next per- 
turbation terms (see the next section about Gaus- 
sian displacements). 



It is important to notice that Eas . 1391 and 1401 are quite 
more complex than the result obtained in a naive way in 
Sec. lIIII of the paper by simply using the continuity equa- 
tion for the local conservation of mass which led basically 
to 

P(k) = fc 2 G(fc). 

We see that respect to this simple approximation, even 
in the small fc limit (i.e., large spatial scales) and finite 
displacements (i.e., finite G(0)), there can be important 
corrections coming both from the granularity of the sys- 
tem, from the internal correlations of the initial mass 
distribution, and from the interplay between these corre- 
lations and those of the displacement field. 
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VIII. CORRELATED GAUSSIAN 
DISPLACEMENT FIELD. 



In this section the effect of a correlated displacement 
field on the correlation properties of a given particle dis- 
tribution is further clarified through the discussion of 
a very important example: the Gaussian displacement 
field. Its importance is twofold: (1) its statistical prop- 
erties are completely determined by the knowledge of 
one and two-point correlation functions (i.e., mean value 
u and correlation matrix G jUI/ (x)); (2) in many appli- 
cations (e.g., initial conditions of cosmological n— body 
simulations) the Gaussian of perturbations (i.e., particle 
displacements) arises as a natural hypothesis basically 
due to arguments based on the central limit theorem |2l| . 
Moreover, in order to clarify better all the concepts in- 
troduced in the previous section, two explicit examples 
of application of a Gaussian displacement field will be 
given: stationary Poissonian and regular lattice initial 
conditions. 

We treat in detail the particular case of a one- 
dimensional spatially stationary point process ni n {x) 
perturbed by a statistically stationary Gaussian dis- 
placement field u{x). The probability density functional 
Q[u{x)] giving the statistical weight of any realization 
u(x) of the field will have the form [25|: 



Q[u(x)] ~ exp 



dxdyu{x)K(x — y)u(y) 



where K(x) = K(—x) is the positive definite correlation 
kernel of the Gaussian displacement field. Without loss of 
generality we have put u = 0. Clearly Eq. 021 is only for- 
mal as, rigorously speaking, the normalization constant 
is infinite in the continuum space and also in the infinite 
volume limit. However we write it to make evident the 
meaning of Gaussianity for a stochastic field. For a more 
general account of Gaussian stochastic fields see |26j, |30| . 

It is simple to show that the PS of the displacement 
field is given by 



G(k) = 



1 



FT[K(x)} 



As we want u(x) to be a well defined continuous stochas- 
tic process, K(x) is to be taken so that the Wiener- 
Khinchin theorem is satisfied j2^], i.e., K(x) is such that 
G(k) > at all k and integrable. The displacement- 
displacement correlation function will be thus given by 



1 



G(x) = u(x Q + x)u(x ) = FT 1 



FT[K(x) 



which is a continuous function if u(x) is a continuous 
stochastic process [2l]. It is possible to show that the 



joint two-displacement PDF can be written as: 

f(u,v;x) 



2tt v /G 2 (0) -G 2 {x) 

G(0)(u 2 + v 2 ) - 2G(x)uv 



(43) 



exp 



2(G 2 (0)-G 2 (x)) 



where G(0) = u 2 = v 2 < +oo. It is also simple to verify 
that, in order to have this joint PDF well defined at all x, 
the correlation function G{x) must satisfy the following 
constraint: 

\G{x)\ < G(0) Vx ^ 0. 

Note that if, as supposed, the Gaussian displacement 
field is a real continuous correlated stationary stochas- 
tic process and, consequently G(x) is a continuous func- 
tion, then lim x _,o G(x) = G(0). This implies then 
f(k, — k; 0) = 1. If instead the displacement field is un- 
correlated (i.e., G(x) = for x ^ but G(0) > 0), 
then, as shown in Sec. EH f(k, -k; 0) = cxp[-fc 2 G(0)] = 
\p(k)\ 2 with p(k) = exp[— fc 2 G(0)/2] being the character- 
istic function of the one-displacement PDF. 

In the case of truly continuous process, by taking the 
FT of Eq. 031 both in u and v, we find: 



f(k, —k; x) = 0(fc; x) = e -k 2 (G(0)-G( x)) 



(44) 



Therefore, by using Eq. 1441 the relation 123 between the 
PS of the particle distribution after the application of 
the displacement field, and its initial correlation function 
£in(x) w iU be 



P(*)= e ~ kG ^ 
2irS(k) . 



dxe~ lkx+k2G ^{l+i ln {x)) 



(45) 



Since at large x both G{x) and S,i n {x) converge to zero, 
the FT in x in Eq. 031 is not well defined and contains a 
Dirac delta function contribution exactly compensating 
the last Dirac delta function term in the same equation. 
This can be made more clear in the following way. The 
Dirac delta function of Eq. 03] can be transformed as: 



2Tv5(k) = 2ne- k2G ^5{k) = e~ k ' G 



(0) 



+ 00 



dxe~ lkx 



Using this relation, Eq.05jcan be rewritten in the follow- 
ing form containing only well defined FTs: 



P(k) 



-k 2 G(0) 



dxe~ lkx (e k2G{x) - l) + 



dxe- lkx+k2Gix ^ m (x) 



(46) 



The small k expansion of this formula can be obtained 
by simply applying to this case Eq. 0U| 
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As already shown above, in d > 1 dimensions the 
more general form of /(u, v;x), even in the limited 
case of spatially stationary Gaussian displacement field 
is more complex. In fact in general not only parallel 
components of u and v can be correlated, but also per- 
pendicular components can be. This leads to have a 
symmetric displacement-displacement correlation matrix 
G A1 i / (x). Once the correlation matrix is given, it is sim- 
ple to show that, for a d— dimensional spatially stationary 
Gaussian displacement field with zero average value, we 
can write: 



0(k;x)=/(k,-k;x) = 

l,d 

exp -^A>)A^)[G^(0)-G MI ,(x)] 



However, in the case G All/ (x) = 5^G(x), the general fea- 
tures of the effect of a Gaussian displacement field can be 
well summarized by the above one-dimensional example. 



A. Gaussian displacement field applied to a 
Poisson particle distribution 

In this section we analyze in detail the effect of a 
Gaussian displacement field as presented above on a one- 
dimensional Poisson particle distribution. The extension 
to higher dimensions is straightforward. 

Let us suppose to have a spatially stationary Poisson 
particle distribution with constant average number den- 
sity tiq > 0. As aforementioned in this case the NCF 
is 



Sjxl 
n 



i.e., hin(x) = and Pi n (k) = By substituting this 
initial NCF into Eq.^ one finds: 

P(k) = — + e- fe2G (°) [°° dxe~ lkx (e~ k2G ^ - l) , 
no J-oo V ' 



which, in the long wave-lengths limit, behaves as 

P(k) ~ — + k 2 \6(k) - fc 2 G(0) 
n L 



(47) 



Since in d = I, G(k) ~ k? with (3 > -1, we see that at 
small k the leading term is l/n = Pi n (k) and the pertur- 
bations introduced by the displacements are at most of 
order k. However, as aforementioned, in d > 2 "critical" 
and dominating perturbations at small k can be intro- 
duced by the displacements. On the other hand, as clear 
from Eg. 1471 we have a perturbation of order k 4 at small k 
due to the finite variance of the displacements G(0) which 
prevails on the perturbation coming from displacement- 
displacement correlations (i.e., from the term containing 
G(k) in Eq.^I) if /3 > 2. We ultimately observe that for 



any d it is not possible to transform a Poisson particle 
distribution into a "superhomogeneous" one by the only 
action of the displacements. In fact, even though cor- 
related, the displacements are stochastic perturbations 
and, consequently, cannot increase the level of large scale 
order of the system. 



B. Shuffled lattice with correlated Gaussian 
displacements 

In this subsection we analyze the effect of the same 
Gaussian displacement field as above on a one dimen- 
sional lattice (i.e., a regular chain of unit mass particles) 
with lattice spacing I (i.e., uq = l/l). Again we focus our 
attention on the large scales (k — > limit). By taking 
the inverse FT of Eq. 1261 the initial NCF in this case is 

€in(x)=l Y, S(x — nl) — l. 

n— — oo 

Therefore Eq. 1451 can be rewritten as 

+00 

P(k) = le- k2G ^ Y e- lknl+k2G ^ - 2ir5{k) . (48) 



The series in Eq. 0H]is not well defined because its argu- 
ment does not converge to zero for |n| — ► oo. By using 
for 27r<5(fc) the chain of identities 



2n5{k) = 2ire- k2G ^ x 



k - 2tt- 



X. . • X 



„ m 

~ T 



— x . • x 



n— — oo 

= le~ k2G ^ x e ~ lknl 



I 



~2ne- kG ^x Y d(k-2Tr 
we can rewrite Eq. 0H| as 



+ oo 



P(k) = le- k G <°) 



^ — iknl 



k 2 G{nl) _ i 



— 00, + OG 

m/0 



(49) 



The first sum is a well defined series and gives in gen- 
eral the smooth part of P(k) converging to 1/no at large 
k. Instead the second one gives the singular delta-like 
Bragg peaks contribution to P(k) due to the initial per- 
fect lattice distribution, each peak being modulated by 
an amplitude e~ k G (°), rapidly decreasing to zero with 
increasing the position k of the peak. In d > 1 it is 
simple to obtain a very similar formula. 

Since this second contribution, due to the initial NCF, 
is exactly zero at all orders in the region —2tt/1 < k < 
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FIG. 5: The figure presents the contribution to Eg. 1491 com- 
ing from the first sum (i.e., excluding the delta-like Bragg 
peaks) to the PS P(k) of a Id "shuffled lattice" obtained 
by applying a correlated Gaussian displacement field with 
G(k) = exp(-|fc|/fc ) with k = 1/4 (hence G(0) = ~ 0.08) 
to a regular chain of particles with unitary lattice spacing 
(1 = 1). Both the theoretical prediction given by Ea. Inland 
the numerical result obtained by a direct simulation with 10 4 
particles are given. The agreement is excellent. Note that 
since in this case Gd(0) > and finite we have that P(k) ~ k 2 
at sufficiently small k. 

2ir/l, the small k behavior of P(k) is completely deter- 
mined by the first sum. This situation is very different 
from the case of an initial Poisson particle configuration. 
In fact, as shown above, in this second case the small k 
properties of the PS are determined mainly by the dis- 
placement field only if it produces "critical" large scale 
correlations (which is possible only in d > 2). Instead, in 
the present case, due to the particular properties of the 
lattice PS, it is always the displacement field to deter- 
mine the large scale correlation properties, i.e., the small 
k PS. We now analyze in detail the small k behavior of 
P(k) in Eg. 1491 with a particular attention to what kind 
of superhomogeneous distributions we can obtain after 
the application of the displacements. By keeping only 
the most important terms of Ea. 1491 at small k, we can 
write: 

P(k) ~ e- k2G ^k 2 G D (k) + y [G {2 \k) - 2G(0)G D (k)} , 

(50) 

plus higher order terms. In Eq. [5D| Grj(fc) is the dis- 
cretized Fourier integral of G(x) with finite integration 
element given by the lattice cell size I, i.e., 

G D (k) = le~ lknl G(nl), (51) 

71= — OO 

and analogously 

G [2] (k) = Yl le- tknl [G(nl)] 2 . 

n— — oo 



Note that, due to the discretization of the Fourier in- 
tegral in a sum in Eq. |^ Go(k) is not exactly equal to 
G(k). This difference depends on how smooth and con- 
stant is G(x) on the length scale I and vanishes when I 
goes to zero. However in general if G(k) ~ k a at small k 
with a < then also Go{k) ~ k a . Instead if a > 0, and 
then G(0) = 0, we can have Gd(0) > (but vanishing as 
k a when / — > 0) as the effect of the discretization of the 
Fourier transform. Another important observation about 
Eq. G>0]is that, as [G(x)\ 2 > at all x, then G^(0) > 0. 
Moreover, as [G(x)] 2 decays faster than G(x) at large x 
then, if Gd(0), is finite also G y D ' (0) is. 

These observations are very important in order to de- 
termine which kind of point processes can be realized 
by perturbing a lattice through a correlated stochastic 
displacement field. By the previous considerations, it is 
straightforward to see that the realization from a lat- 
tice of a stochastic particle distribution with P(k) ~ k 13 
at small k with /3<2(l</3<2inrf=l, and 
2 — d < (3 < 2 va general d dimensions) is a very sim- 
ple task: it is enough to take a Gaussian displacement 
field such that G(k) ~ k a with a = (3- 2 < (see Figs.0 
andHjJ). Instead, obtaining a P(k) with (3 > 2 is a more 
difficult task. First of all one has to take a G(x) such that 
Gd(0) = 0. It is possible to see that this requirement is 
a sort of stochastic expression of the total momentum 
conservation law in the system, that is, a conservation 
"in average" |29j. Moreover in order to have 2 < (3 < 4 
one has to have Gr>(k) ~ k 13 ^ 2 at small k. If instead 
Gn(k) ~ k a with a > 2 then (3 — 4 is obtained in all 
cases. The case (3 > 4 is not permitted. This means that, 
by perturbing a lattice with a (Gaussian) stochastic dis- 
placement field, the most "superhomogeneous" realizable 
point process has a PS P(k) ~ k 4 at small k, e.g., k e is 
forbidden. This is a strong limitation. In fact, as afore- 
mentioned, from one side the PS of the lattice can be 
considered as behaving ~ k°° around k = (it is com- 
pletely flat). From the other side we have just shown 
that an exponent (3 > 4 of P(k) is never realizable when 
one applies to any initial particle distribution a stochas- 
tic displacement field. This means that independently of 
the kind of the displacement field, as long as a smooth 
G(k) around k = is considered, there is a "minimal 
level of disorder" injected in the system measured by the 
exponent 4 in the small k behavior of P(k). This shows 
how difficult is for example to build a spatially stationary 
stochastic point process such that P(k) ~ k 6 at small k. 



It is possible to extend all these results to higher di- 
mension and to non-Gaussian displacement field, but we 
think that the example of the one-dimensional Gaussian 
displacement field is enough to enlighten the effects and 
the limitations of a general stochastic and correlated dis- 
placement field on a given point process statistically in- 
dependent of it. 
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FIG. 6: Power spectrum P(k) (including also the Bragg peaks 
contribution) of a shuffled lattice similar to the previous fig- 
ure but now with G(k) = Aexp(— k/k c )/(k + a) with k c = 5, 
a = 1CT 4 , and A = 1/20 (hence G(0) = u? = 0.13). The 
agreement between the theoretical prediction and the numer- 
ical result is again excellent. Note that now, as G(k) ~ k^ 1 
at small k (but larger than a), the PS P(k) ~ fe in the same 
range. 
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FIG. 7: The figure presents a comparison between the PSs of 
(i) a Poisson point process, (ii) a shuffled lattice with uncor- 
rected displacements with a finite variance, and (iii) a shuf- 
fled lattice obtained by applying a correlated displacement 
field and showing a "critical" behavior at sufficiently small 
k. In particular all the particle distributions are in d = 1 
and have the same mean density no (i.e., the same specific 
volume 1/no). In the case (ii) the finite variance u 2 < +oo 
implies P(k) ~ k 2 at small k (see Eg. 1231 . Instead in (iii) the 
applied correlated displacement field for a k 2ir/na has 
G{k) ~ k~ 4 , and then, from Eq. EH -P(fc) ~ k 2 . For fc < a 
(which is below the minimal k of the figure) G(k) is cutoff 
to a constant value, in fact in any d the PS of any proper 
stochastic process has to satisfy lim fc ^ k d P(k) — 0. 



IX. CONCLUSIONS 

In this paper we have discussed the effect of a stochas- 
tic displacement field on a spatial distribution of point- 
particles with identical mass (i.e., point process) in the 
hypothesis of statistical independence between the two. 
In particular we have studied rigorously the changes in- 
duced in the two-point correlation properties of the par- 
ticle distribution by the displacements. In this way we 
have seen that usual naive approaches to this problem 
leads to approximations (which may not be appropriate 
in some important cases) neglecting the contributions 
coming from either the initial correlation properties of 
the particle distribution or the failure of the small dis- 
placements approximation. 

We have distinguished the two cases of displacement 
fields with and without spatial correlations, giving for 
both cases the rigorous equations of transformation of 
the PS of the particle distribution. 

Our main interest concerns a detailed study of the kind 
of large scale correlations that the displacement field can 
"inject" into the system. In particular we are interested 
to know if long range positive spatial correlations in the 
system can be obtained by the application of a suitable 
displacement field to an initially short range correlated 
particle distribution. In this respect we have found that 
in d > 3 it is possible to obtain particle distributions with 
long range correlations (i.e., with a covariance function 
with a positively diverging integral) by applying for in- 
stance either to a completely uncorrelated homogeneous 
Poisson point process or to a regular lattice of particles a 
Gaussian displacement field with sufficiently long range 
displacement-displacement correlations independently on 
the variance of the single displacement. This is a very im- 
portant point, in fact we can think to perturb a lattice 
with such a displacement field with a mean squared dis- 
placement much smaller than the lattice spacing, and, 
independently of this, to obtain strong large scale typical 
fluctuations AN(R) of the number of points N(R) con- 
tained in a randomly placed sphere of radius R, growing 
as AN(R) ~ R a with a > d/2, whereas a = d/2 for the 
Poisson distribution and a = (d — l)/2 for the regular 
lattice. Such particle distributions would look very simi- 
lar respectively to the original Poisson or lattice particle 
distribution at the small scales (i.e., locally), but would 
show much larger (and more rapidly growing with the 
distance) fluctuations beyond a sufficiently large scale. 

Another related problem is to study how the long range 
order of a regular lattice array of particles is perturbed 
by the action of the applied displacement field. In partic- 
ular we have studied between which limits the perturbed 
lattice stays in the so-called superhomogeneous class of 
point processes. That is, we have found the limits in 
which the sub-Poissonian character of the mass (i.e. the 
particles number) fluctuations on sufficiently large scales 
is conserved. Anyway we have also seen that any truly 
stochastic displacement field always injects into the parti- 
cle system at least a minimal degree of disorder measured 
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by the a = 4 of the exponent of the final PS of the system 
P(k) ~ k a at small k. 

We think that all these results can be of very practical 
importance in many physical and more largely scientific 
applications. For instance, this is the case of the so called 
n-body cosmological simulations. In this case the super- 
position of a suitable stochastic displacement field to a 
pre-initial particle distribution (e.g., a lattice of identical 
particles) is the usual way to prepare the initial condi- 
tions for numerical simulations to study the problem of 
the dynamic gravitational clustering of the matter lead- 
ing to the formation of large scale structures (e.g., galax- 
ies) . These so built initial conditions should represent the 
spectrum of small primordial mass fluctuations predicted 



by theoretical models (e.g., CDM models). 
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In the case of ergodicity it can be taken also to be a 

volume average in the infinite volume limit. 

E.g., the case in which p(u) — S(u — uo)/2 + <5(u + uo)/2 

for which p(k) = cos(k • uo) and then |p(k)| = 1 for all k 

such that k • uq = nir with n any integer. 

As |p(0)| = 1, for k — instead we find singularly 

5 m (0) = Sin(0) at any m. This creates in general an 

asymptotic isolated discontinuity at k — 0. However this 

does not matter for the real space correlation properties 

as the Fourier transform of the PS is insensible to that. 



